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The 0(a 2 ) massive operator matrix elements for unpolarized and polarized heavy flavor 
production at asymptotic values Q 2 >> m 2 are calculated in Mellin space without 
applying the integration-by-parts method. We confirm previous results given in Refs. [4, 
5], however, obtain much more compact representations. 

1 Introduction 

The heavy-flavor corrections to deeply inelastic structure functions are very important for 
the range of small values of x and do contribute there on the level of 20-40%. They have to 
be known at the same level of accuracy as the light-flavor contributions for precision mea- 
surements of Aqcd [1] and the parton distributions. The next-to-leading order corrections 
were given semi-analytically in [2] for the general kinematic range. Fast and accurate imple- 
mentations of these corrections in Mellin-spacc were given in [3]. In the region Q 2 >> m 2 , 
the heavy flavor Wilson coefficients were derived analytically to 0(a 2 ) [4,5]. Here Q 2 de- 
notes the virtuality of the gauge boson exchanged in deeply-inelastic scattering and m is the 
mass of the heavy quark. In this note we summarize the results of a first re-calculation of 
the operator matrix elements (OMEs) in [6,7]. The calculation is being performed in Mellin- 
space using harmonic sums [8, 9] without applying the integration- by-parts technique. In 
this way, we can significantly compactify both, the intermediary and final results. We agree 
with the results in [4,5]. The unpolarized and polarized 0(a 2 ) massive OMEs can be used 
to calculate the asymptotic heavy-flavor Wilson coefficients for F 2 (x,Q 2 ) and gi(x,Q 2 ) to 
0(a 2 ) [4-7], and for F L (x,Q 2 ) to 0(a 3 s ) [10]. 

2 The Method 

In the limit Q 2 >> m 2 the heavy quark contributions to the twist-2 Wilson coefficients are 
determined by universal massive operator matrix elements (i\Ai\j) between partonic states. 
The process dependence is due to the corresponding massless Wilson coefficients [11]. This 
separation is obtained by applying the renormalization group equation(s) to the (differential) 
scattering cross sections, cf. [4]. In this way all logarithmic and the constant contribution 
in m 2 IQ 2 can be determined. The operator matrix elements are calculated applying the 
operator insertions due to the light-cone expansion in the respective amplitudes. One obtains 
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with (g> denoting the Mellin convolution. The OMEs contain ultraviolet and collinear di- 
vergences. The collinear singularities are absorbed into the parton distribution functions 
while the ultraviolet divergences are removed through rcnormalization. To 2-loop order, the 
renormalized OMEs read : 
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and similar for the quarkonic contributions. Here, fi 2 denotes the factorization and rcnormal- 

(k—l) 

ization scale, are the feth loop splitting functions and /3q denotes the lowest expansion 

coefficient of the /3-function. and arc the O(e ) rcsp. 0(e)-terms in the expansion 
of the OME, which form the main objective of the present calculation. 

3 Results 

We calculated the massive operator matrix elements both, for the gluon-heavy quark and 
light-heavy quark transitions in the flavor non-singlct and singlet cases, for unpolarized and 
polarized nuclcon targets. 

The constant contribution to the unpolarized and polarized OMEs for the transition 
g — » Q are : 
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Here Pi{N) denote polynomials given in [6,7]. The corresponding quarkonic expressions 
are given in [6,7]. The integrals were performed using Mcllin-Barnes techniques [12,13] 
and applying representations in terms of generalized hypergcomctric functions. The results 
were further simplified using algebraic relations between harmonic sums [14]. Furthermore, 
structural relations for harmonic sums [15], which include half-integer relations and differ- 
entiation w.r.t. the Mellin variable N, lead to the observation that the OMEs above depend 
only on two basic harmonic sums : 

Si(N), S- 2 ,i(N) ■ 

We expressed S-2,i(N) in terms of the Mellin transform M[Li 2(2) / (1 + x)](N) in the above. 
Here f3(N) = (1/2) • [ip((N + l)/2) — ip(N/2)]. Previous analyzes of various other space- 
and time-like 2-loop Wilson coefficients and anomalous dimensions including also the soft 
and virtual corrections to Bhabha-scattering [15a, 16], showed that six basic functions are 
needed in general to express these quantities : 

Si(N), S ±2 ,i(N), S_3,i(JV), S±2,i A (N) . 

None of the harmonic sums occurring contains an index { — 1} as observed in all other cases 
being analyzed. 

Comparing to the results obtained in Refs. [4, 5] in x-space, there 48 functions were 
needed to express the final result in the unpolarized case and 24 functions in the polarized 
case. 

To obtain expressions for the heavy flavor contributions to the structure functions in 
a; -space, analytic continuations have to be performed toJVeC for the basic functions given 
above, see [15, 17, 18]. Finally a (numeric) contour integral has to be performed around the 
singularities present. 



4 Conclusions 

We calculated the unpolarized and polarized massive operator matrix elements to 0(a 2 ), 
which are needed to express the heavy flavor Wilson coefficients contributing to the deep- 
inelastic structure functions F2,gi and Fl to 0(a 2 ) resp. 0(a 3 s ) in the region Q 2 » m 2 . 
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The calculation was performed in Mellin space without using the integration-by-parts tech- 
nique, leading to nested harmonic sums. We both applied representations through Mcllin- 
Barnes integrals and generalized hypergeometric functions. In course of the calculations, a 
series of new infinite sums over products of harmonic sums weighted by related functions 
were evaluated, cf. [6,7]. These representations were essential to keep the complexity of 
the intermediary and final results as low as possible. Furthermore, we applied a series of 
mathcmatic relations for the harmonic sums to compactify the results further. We confirm 
the results obtained earlier in Refs. [4, 5] by other technologies. 
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